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Abstract 

Background: Plant color variation is due not only to the global pigment concentration but also to the proportion 
of different types of pigment. Variation in the color spectrum may arise from secondary modifications, such as 
hydroxylation and methylation, affecting the chromatic properties of pigments. In grapes {Vitis vinifero L), the level 
of methylation modifies the stability and reactivity of anthocyanin, which directly influence the color of the berry. 
Anthocyanin methylation, as a complex trait, is controlled by multiple molecular factors likely to involve multiple 
regulatory steps. 

Results: In a Syrah x Grenache progeny, two QTLs were detected for variation in level of anthocyanin methylation. 
The first one, explaining up to 27% of variance, colocalized with a cluster of Myb-type transcription factor genes. 
The second one, explaining up to 20% of variance, colocalized with a cluster of 0-methyltransferase coding genes 
(AOMT). In a collection of 32 unrelated cultivars, MybA and /\0M7 expression profiles correlated with the level of 
methylated anthocyanin. In addition, the newly characterized A0MT2 gene presented two SNPs associated with 
methylation level. These mutations, probably leading to a structural change of the A01\/1T2 protein significantly 
affected the enzyme specific catalytic efficiency for the 3'-0-methylation of delphinidin 3-glucoside. 

Conclusion: We demonstrated that variation in methylated anthocyanin accumulation is susceptible to involve 
both transcriptional regulation and structural variation. We report here the identification of novel AOMT variants 
likely to cause methylated anthocyanin variation. The integration of QTL mapping and molecular approaches 
enabled a better understanding of how variation in gene expression and catalytic efficiency of the resulting 
enzyme may influence the grape anthocyanin profile. 



Background 

Anthocyanins represent a major group of the flavonoid 
family, consisting of small water-soluble molecules stored 
in vacuoles. As widespread plant secondary metabolites, 
they are responsible for red and blue colors of many 
plant tissues [1]. Anthocyanins, as major pigments, play 
important roles in plant reproduction by attracting polli- 
nators and seed dispersers and in protection from photo- 
oxidative stress [2]. In crops, ornamental plants, and 
fruits, they constitute a key trait for the aesthetic quality 
of the product. In grape, they play a crucial role as they 
participate in both wine color and organoleptic proper- 
ties due to their complex interactions with other phenolic 
compounds, as well as with proteins and polysaccharides 
[3]. 



* Correspondence: ageorges(a)supagro.inra.fr 

^INRA, UMR 1083 Sciences pour I'OEnologie, 34060 Montpellier, France 
Full list of author information is available at the end of the article 



The anthocyanins identified in Vitis vinifera are 3-0- 
monoglucosides (3-glc) and 3-0-acyl monoglucosides 
derived from the five main anthocyanidins-delphinidin 
(Dp), cyanidin (Cy), peonidin (Pn), petunidin (Pt) and 
malvidin (Mv)-, which differ from each other in the 
number and position of the hydroxyl and methoxyl 
groups located on their B-ring (Figure 1): delphinidin and 
cyanidin being unmethylated, peonidin and petunidin 
mono-methylated and malvidin di-methylated. Methyla- 
tion by S-adenosyl-L-methionine (SAM) dependent-0- 
methyltransferase plays an important role in modifying 
the structure of plant secondary metabolites. For antho- 
cyanins, methylation stabilizes the phenolic B ring, thus 
reducing the reactivity of the overall molecule [4], and 
increases water solubility, thus reinforcing its color prop- 
erties [5,6]. The color of anthocyanins greatly depends on 
the number of hydroxyl groups: the larger the number of 
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Figure 1 Chemical structure of anthocyanins from grape berries {Vitis vinifera L). 
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groups, the bluer the color. On the other hand, the O- 
methylation of anthocyanins has a reddening effect. 

The anthocyanin biosynthesis pathway has been exten- 
sively characterized in numerous plant species [7,8]. In 
grapevine, the core structural genes of the anthocyanin 
pathway leading to the synthesis of Dp 3-glc and Cy 3-glc 
have been cloned and characterized [9,10]. Recently, two 
anthocyanin 0-methyltransferases (AOMT and FAOMT) 
able to methylate anthocyanins in vitro and in vivo have 
been identified [11,12]. Anthocyanin biosynthesis is con- 
trolled by Myb-bHLH transcriptional complexes that acti- 
vate the expression of structural genes involved in the late 
steps of the pathway [7,13-16]. In particular, the VvMybA 
genes induce the transcription of UFGT and AOMT in 
colored tissues [17-19] and were shown to control antho- 
cyanin accumulation in grape [20,21]. While anthocyanin 
methylation in grapes is altered by environmental and cul- 
tural conditions [22,23], it is also strongly affected by 
genetic factors. Indeed, both the total amount of anthocya- 
nins and the relative abundance of single anthocyanins are 
extremely variable among red-to blue-skinned cultivars. 
In spite of this variability, mono- and di- methyl derivatives 
are largely predominant [24,25], malvidin 3-glucoside 
being the main anthocyanin in most cultivars and peoni- 
din 3-glucoside usually being less abundant [26,27]. Genes 
encoding flavonoid 3'-and 3'5'-hydroxylases involved in 
anthocyanins hydroxylation [28,29] or genes encoding O- 
methyltransferase involved in anthocyanins methylation 
[11,12,30] can potentially generate differences in anthocya- 
nin composition. The availability of the full genome 
sequence [31,32] as well as broad genetic resources will 
help us to further clarify the mechanisms involved in 
anthocyanin biosynthesis in grape. 



The objective of this study was to understand the 
architecture of the genetic control of anthocyanin 
methylation in grape and to improve our understanding 
of the complex factors involved in the regulation of 
anthocyanin methylation. Two QTLs for anthocyanin 
methylation variation were identified: one colocalized 
with a MybA gene cluster and the second colocalized 
with a cluster of three putative AOMT genes. This work 
aims to validate the role of the latter in anthocyanin 
methylation and to identify the molecular factors under- 
lying genetic variation. The combination of forward 
genetic approaches (pedigree and population-based) 
with functional genomics and enzymology allowed the 
description of the complex regulation of an anthocya- 
nin- 0-methyltansferase gene cluster as a major player 
for the control of anthocyanin methylation in grape. 

Results 

Variation in methylated anthocyanin content 

In order to investigate the genetic bases of anthocyanin 
methylation, three plant collections were analyzed: (1) a 
Syrah x Grenache (SxG) progeny of 191 genotypes repli- 
cated in two blocks, (2) a collection of 50 colored geno- 
types from the Vassal germplasm collection (Coll 50 cv, 
INRA, Domaine de Vassal, France), and (3) a subset of 32 
genotypes selected to optimally represent the diversity in 
anthocyanin methylation level of the original 50 indivi- 
duals sample with a smaller number of individuals (Core- 
Coll 32 cv). Anthocyanin composition differed dramati- 
cally both in the SxG progeny (Additional file 1) and in 
the collection of natural diversity (data presented for 
Core-Coll 32 cv. Additional file 2). After a /o^- transfor- 
mation, distribution of methylated anthocyanin (MT) 
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was normal in the three samples (P-val > 0.05). The SxG 
progeny showed extensive variation in levels of MT, ran- 
ging from 50% to 97% of methylated anthocyanins with a 
MT variance of cj^mt SxG = 0-42 (Figure 2). The collec- 
tion of 50 colored grape varieties (Coll 50 cv) showed 
variation ranging from 31% to 97% of methylated antho- 
cyanins with a MT variance of (T^mt so cv = 0-58. The var- 
iance was increased to cj^mt 32 cv = 0-89 in the Core-Coll 
32 cv. 

QTL mapping of anthocyanin methylation level 

QTL detection was performed either on the two field 
blocks taken separately or combined using a multi-trait 
procedure. In both cases, this allowed the consistent iden- 
tification of two QTLs for the anthocyanin methylation 
level. The first QTL, explaining 27% of the phenotypic 



variance for anthocyanin methylation, was located on 
Linkage Group 2 (LG2) between the markers VMC5G7 
and VMC8C2 (Figure 3). This QTL was detected on the 
Syrah parental map and supported with a maximum LOD 
of 7.49 for the multi-trait detection. It was also detected 
on both blocks taken separately. This QTL, which coloca- 
lized with a QTL previously identified as controlling the 
content of anthocyanins in grape berry [20], corresponded 
to the VvMybA genes cluster, thus controlling together 
anthocyanin methylation and total accumulation. The sec- 
ond QTL was identified on LGl of the Grenache parental 
map between markers VMC9E3 and VMC8E8, explaining 
up to 20% of the phenotypic variance with a maximum 
LOD of 6.96 (Figure 3). In the present study, we focused 
on this second QTL, which was specific to anthocyanin 
methylation. The marker pOMT2, the closest to the LGl 
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Figure 2 Distribution of the MT variable in the samples used in QTL and association mapping. The /W7 variable is expressed as the 
logarithm of the ratio of methylated/non-methylated anthocyanins in mg of anthocyanin per g of fresh berry skin. SxG: 191 individuals from the 
Syrah x Grenache progeny used in QTL mapping. Coll 50 cv: random sampling of 50 colorful grape cultivars from the original 
agromorphological core-collection (Barnaud et al. [33]). Core-Coll 32 cv: core-collection designed from the Coll 50 cv to optimize the capture of 
A/IT variability. 
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Figure 3 Position and linkage intensity of the QTL associated with the MT trait. The two QTLs were detected on the SXG mapping 
population on the Grenache parental map and on the Syrah parental map, respectively. Confidence intervals correspond to the 90% and 99% 
probability of presence of the QTLs after 1000 bootstrap resamplings. The QTLs presented were found for the A and B blocks and using multi- 
trait detection. 



QTL, showed a difference of 11.8% between homozygous 
combination of alleles leading to high level of MT and het- 
erozygous combination at pOMT2 leading to low level of 
MT (Additional file 3). 

Candidate genes selection in the QTL on LG1 

The physical interval of 1 Mbp surrounding the maxi- 
mum LOD of the QTL of LGl encompassed 46 predicted 
unigenes with an identity score superior to 5 (Additional 
file 4). Among this unigenes set, a cluster of three genes 
encoding putative anthocyanin 0-methyltransferases 
(AOMT) within a 33 Kbp region was identified, referred 
to hereafter as VvAOMTl, VvA0MT2 and VvAOMT?> 
(Figure 4). VvAOMTl encoded a protein that is identical 
to the AOMT from Syrah shown to methylate anthocya- 
nin both in vitro and in planta [11] and similar to the 



FAOMT from Cabernet Sauvignon shown to methylate 
anthocyanin in vitro [12]; therefore, we narrowed the 
study to these three candidate genes. The three genes 
have similar structures, with five small exons ranging 
from 51 bp to 294 bp and four introns (Figure 4). The 
deduced amino acid sequences of the three OMTs 
showed high sequence similarity (Additional file 5) with a 
minimum pairwise identity score of 90% at the genomic 
level and of 96% at the transcript level. In addition, there 
was no other protein showing more than 20 amino acids 
with 80% similarity to AOMT in the grape genome, 
showing evidence for VvAOMT being a single family of 
paralogs. The corresponding transcripts were 708 bp 
long and very similar to VvAOMTl, All three genes were 
thus potentially functional. The 3' non-coding ends of 
VvAOMr genes appeared less conserved than the coding 



Fournier-Level et al. BMC Plant Biology 201 1, 11:179 
http://www.biomedcentral.eom/1471-2229/1 1/179 



Page 5 of 14 



VvAOMTS VvA0MT2 
GSVIVT01 01 0466001 GSVIVT01 01 0468001 

8.3Kbp 



VvAOMn 
GSVIVT01 01 0469001 



19.6Kbp 



I I 
I I 



I I 
I I 



I I 
I I 



i 1.45Kbp 1 
^ ►! 




2.17Kbp 
\< ^ 



2.14Kbp ] 
< ►! 

Eh-Chd 



850bp Sequenced on genomic DNA for the Col 50cv 



700bp Sequenced on coding DNA for the CoreCol 32cv 



Figure 4 Genomic structure of thie open reading frames of the \/v>AO/Wr genes cluster on chromosome 1. The continuous lines represent 
sections drawn at tine riglit scale; dashed lines represent change in scale. Lines between dots represent the sequenced fragments. The filled 
boxes indicate the exons. Data correspond to the 12X genome annotation available from the Genoscope Grape Genome Browser database 
http://www.genoscope.cns.fr/externe/GenomeBrowser/Vitis/. 



region, and the only robust distinction between the three 
isogenes thus relied on divergent non-coding 3' ends. Iso- 
gene-specific primers could be designed for VvA0MT2 
and VvA0MT3 cDNA, but not for VvAOMTl cDNA, 
due to sequence homology with the two other isogenes 
and a shorter 3' end. Using specific primers for 
VvAOMTS, no expression could be detected in mature 
berries of Syrah, Joubertin and Petit Bouschet, suggesting 
that this isogene was not expressed in mature red berries 
(Additional file 6). Although this does not demonstrate 
the non-functionality of VvAOMTS, it justified the 
further consideration of VvAOMTl and VvA0MT2, as 
candidate genes in association genetic analysis both 
expressed in mature berries. 

Effect of the VvAOMT and VvMybA expression on 
methylated anthocyanins 

To assess the relative influence of VvAOMT and 
VvMybA genes, which potentially underlie the QTLs on 
LGl and LG2, respectively, we measured their respective 
transcript accumulation in the Core-Coll 32 cv, maximiz- 
ing the diversity for methylated anthocyanins. We ana- 
lyzed the level of expression of VvAOMT (eAOMT), 
which combines the expression of the three isogenes 
VvAOMTl, VvAOMTl and VvA0MT3 and of VvMybA 
(eMybA), which combines the expression of the two iso- 
genes VvMybAl and VvMybA2. Both eAOMT and 
eMybA appeared to follow continuous variation among 
the 32 cultivars, showing no evidence of regulation by a 
single major polymorphism. eAOMT appeared to be 



highly correlated with the methylated anthocyanin con- 
tent of the berry skin (p = 0.58, P-val < 0.001; Table 1). 
eMybA appeared correlated with total anthocyanin con- 
tent, but not with the methylated anthocyanin content 
(Table 1). Finally, the methylated anthocyanin content 
was better predicted when the effects of eAOMT and 
eMybA were considered jointly in a multivariate regres- 
sion model (Table 2), suggesting that the expression level 
of both genes had an effect on anthocyanin methylation, 
probably with a direct influence of eAOMT 2ind an indir- 
ect influence of eMybA, 



Table 1 Correlation between the anthocyanin content 
and VvMybA, VvAOMT and specific VvA0MT2 expression 





Total Antho 
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MT 


eAOMT eA0MT2 


Total Antho 


1 
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0.48 


1 








0.042* 








MT 


-0.01 


-0.17 


1 






0.9636 


0.35 






eAOMT 


0.31 


0.27 


0.58 


1 




0.0812 


0.1348 


0.0005** 




eA0MT2 


0.04 


-0.04 


0.02 


0.40 1 




0.8356 


0.8182 


0.93 


0.0229* 



Total Antho: total anthocyanin content, MT: level of methylated anthocyanin, 
eMybA: global VvMybA expression level, eAOMT: global i/v/AOMT expression 
and eA0MT2: specific VvA0MT2 expression In the Core-Coll 32 cv sample. The 
Pearson's p values are presented In plain text while their associated P-values 
are shown In Italic {* < 5% and ** < 1%) 
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Table 2 Effects of VvAOMT and VvMybA global 
expression on the MT variable in the Core-Coll 32 cv 



N of variables ajusted 


BIC 


Variables in the model 




1 0.31 


-23.17 


eOMT 




2 0.41 


-26.72 


eOMT 


eMybA 






0.31 


0.1 



Both variables [eAOMT and eMybA, respectively) show independent significant 
effect in a multiple regression test, the number in italic corresponds to the 
partial R-squared of each effect tested 



Association mapping with VvAOMTI and VvA0MT2 

To isolate the particular effect of the two VvAOMT c^Ln- 
didate genes, we partially sequenced both genes in the 
collection of 50 colored varieties (genomic DNA frag- 
ments corresponding to 700 bp of 3'end of VvAOMTI 
and 850 bp of 3'end of VvA0MT2; Figure 4). This 
allowed detecting polymorphisms on genomic DNA in 
the two last exons and in the 3'UTR specifically for both 
VvAOMT candidates (high sequence homology and 
shorter 3'UTR for VvAOMTI did not allow specific 
amplification at the complete coding sequences). On 
VvAOMTI, we identified 15 SNPs and 2 InDels, and, on 
VvA0MT2, we identified 19 SNPs and 1 InDel, giving 
rise to 9 non-synonymous mutations. We tested the asso- 
ciation of these polymorphisms with the level of methy- 
lated anthocyanins using the EMMA software [34]. Six 
SNPs were identified as associated with variation in 
anthocyanin methylation in VvA0MT2 and none in 
VvAOMTI, The significantly associated polymorphisms 
with P-value < 0.05 and the test statistics are presented 
in Table 3. SNP Y522 did not change the amino acid 
sequence of the gene, and the five other associated poly- 
morphisms were located in the 3'UTR part of the 
VvAOMTI gene, thus not affecting the protein structure. 
All six associated SNPs showed significant LD, with one 
LD block, including SNPs Y522, Y717, Y741 and W774 
(pairwise r ^ = 0.5-1), and another LD block, including 
M843 and InDel909 (pairwise r ^= 0.6). However, given 
the high level of polymorphism at both loci, we expected 
that if there was functional variation elsewhere in the 



VvAOMTI sequence, we would have genotyped at least 
one SNP significantly associated through LD with the 
causal variation. In addition to the two blocks of LD 
reported for VvA0MT2, two blocks of high LD were 
identified on VvAOMTI (Additional file 7). However, no 
intergenic LD was significant between VvAOMTI and 
VvA0MT2, showing a rapid decay of LD in the genomic 
fragment including these genes. As VvA0MT2 was the 
only isogene showing polymorphism association with 
anthocyanin methylation level, we only retained 
VvA0MT2 as candidate gene. 

Specific expression of VvA0MT2 

Isogene-specific primers were designed to measure the 
specific expression of the candidate gene VvAOMTI in 
(1) MybAl transformed grapevine hairy roots accumu- 
lating anthocyanin [18] and in (2) the Core-Coll 32 cv 
(Additional file 2). As we did not succeed in designing 
specific primers for AOMTl, we could only compare the 
global expression eAOMT {VvAOMTI plus VvAOMTI) 
using non specific primers with the specific expression 
VvAOMTI, In hairy roots over-expressing MybAl, the 
expression of VvA0MT2 showed a 5.5-fold increase 
compared to the control, while the level of eAOMT, 
showed a 75-fold increase (Additional file 8). Although 
the effect of VvMybAl on VvA0MT2 cannot be 
neglected, these results suggest that the effect of 
VvMybAl on VvA0MT2 is Umited. In the 32 cultivars 
subcore-coUection, no significant correlation was found 
between the expression level of VvA0MT2 and the ratio 
of methylated anthocyanins. Moreover, no significant 
effect of eMybA was observed on the expression of 
VvA0MT2 (Table 1). However, even if the specific 
expression level of VvA0MT2 is correlated with global 
eAOMT (p = 0.40, P-val = 0.023) and has a significant 
regulatory effect of eAOMT on the level of methylation, 
in the absence of correlation between VvA0MT2 
expression and MT, we concluded that VvA0MT2 
expression had no direct regulatory effect on anthocya- 
nin methylation. 



Table 3 Association genetics tests performed on VvAOMTI polymorphisms for the level of methylated anthocyanin 

Samples Marker df Marker F test P-value Model Marker 

Coll 50 cv Y522 1 12.51 0.0041 0.61 0.41 

Y717 1 5.67 0.0347 0.46 0.26 

Y741 1 10.21 0.0085 0.51 0.46 

W774 1 5.67 0.0347 0.46 0.26 

M843 2 5.75 0.0195 0.61 0.41 

lnDel909 2 4.22 0.0435 0.55 0.35 

Core-Coll Y127 2 7.80 0.0039 0.53 0.43 

32 cv W388 2 7.80 0.0039 0.53 0.43 

Y741 2 7.80 0.0039 0.53 0.43 



Tests were performed using genomic DNA sequence fragments in the Coll 50 cv and coding DNA sequence in the Core-Coll 32 cv. Only the significant tests are 
shown. Df: degree of freedom, R^: percentage of explained variance 
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Sequence analysis of the VvA0MT2 gene 

We fully sequenced the VvA0MT2 cDNA of the Core- 
Coll 32 cv and identified 30 polymorphisms in the ORF, 
including 26 SNPs and 4 InDels (Additional file 9). As the 
genomic and cDNA fragments sequenced for VvA0MT2 
partially overlapped, 8 were common to those already 
identified in the 50 cultivars sample, including the 4 most 
associated SNPs presented in Table 3. This new poly- 
morphism set was tested again for association with methy- 
lated anthocyanin level in the Core-Coll 32 cv and only 
three associations were found to be significant: Y741 was 
again strongly associated, as well as two new SNPs, Y127 
and W388. These three SNPs showed complete Unkage 
disequilibrium (r ^ = 1) and their effects were thus con- 
founded. In addition, the VvOMT2 sequence for the par- 
ents of the SxG population showed to be homozygous for 
the strong allele of W388 and Y741 in the Syrah parent 
and heterozygous in the Grenache parent, following the 
same segregation as the pOMT2 marker. This supported 
that the Grenache QTL on LGl could be attributed to the 
segregation of VvOMT2 alleles in the mapping population. 

Finally, in the best model including only expression pro- 
files, MOMr explained 31% of the variance while the sin- 
gle SNP Y127 or W388 in the association model explained 
43% of variance for the same sample (Table 2 and Table 
3). This strong effect was associated with a highly signifi- 
cant difference in anthocyanin methylation between the 
groups carrying the strong or the weak allele at SNPs 
Y127 and W388 (MT being on average 94.7% and 76.5% 
respectively). The SNPs Y127 and W388 corresponded to 
non-synonymous changes affecting the primary structure 
of the AOMT protein. Y127 led to a substitution from H 
to Y in position 43 of the amino acid chain and W388 led 
to a substitution from S to T in position 130 of the amino 
acid chain. Y127 appeared to generate a mutation within 
the second a-helix of the predicted protein secondary 
structure and W388, within the sixth predicted a-helix. 

Characterization of the allelic variants of VvA0MT2 

Two grape varieties were selected for the homozygous pre- 
sence of contrasting SNPs Y127 and W388 at the 
VvA0MT2 locus, representing two allelic variants of the 
corresponding AOMT2 protein. 'Petit Bouschet x Aramon 
n°4' (PB) presented the weak alleles T at Y127 and A at 
W388 and 'Papadiko' (PK) presented the strong alleles C 
at Y127 and T at W388. No other difference was detected 
in the coding region. Enzymatic tests performed with PB- 
AOMT2 and PK-AOMT2 indicated that both enzymes 
catalyzed in vitro the 3' 0-methylation of anthocyanins 
with a catechol B-ring (3',4' di-OH) and the 3' and 5' O- 
methylation for those showing a pyrogallol (3',4',5' tri-OH) 
B-ring. However detailed analysis of PB-AOMT2 and PK- 
AOMT2 activity revealed substrate-specific differences in 
catalytic efficiency (Table 4). Indeed, when Cy 3-glc was 



Table 4 Substrate specificity of the recombinant A0MT2 
proteins from Papadiico (PK-A0MT2) and from Petit- 
Bouschet X Aramon n°4 (PB-A0MT2) 



Enzyme Substrate Products relative amounts (%) 



PK-A0MT2 


Dp 3-glc 


Dp 3-glc 


51.1 (± 6.9) * 






Pt 3-glc 


40.8 (± 4.0) * 






Mv 3-glc 


8.1 (± 3.3) 




Cy 3-glc 


Cy 3-glc 


58.4 (± 12.2) 






Pn 3-glc 


44.9 (± 15.4) 


PB-A0MT2 


Dp 3-glc 


Dp 3-glc 


63.6 (± 11.1) * 






Pt 3-glc 


31.2 (± 8.6) * 






Mv 3-glc 


5.3 (± 2.6) 




Cy 3-glc 


Cy 3-glc 


60.6 (± 9.3) 






Pn 3-glc 


39.4 (± 9.2) 



Reaction were essayed in vitro and products were identified according to their 
UV visible absorption spectra and retention time. Quantifications were based 
on peak areas at 520 nm. The data presented are the mean of 10 experiments 
performed with 3 independent protein preparations (± standard deviation), * 
indicates significance in ANOVA test (P-values < 0.05). Dp, delphinidin; Cy, 
cyanidin; Pt, petunidin; Pn, peonidin; Mv, malvidin; glc, glucoside 

used as a substrate, both PB-AOMT2 and PK-AOMT2 
appeared to catalyze equally well the methylation of Cy 3- 
glc into Pn 3-glc. By contrast, when Dp 3-glc was used as 
a substrate, PB-AOMT2 showed a significantly reduced 
efficiency compared to PK-AOMT2 in producing Pt 3-glc 
(mean difference in relative final Pt 3-glc proportion of 
9.7%, P-val = 0.02, Figure 5). Finally, both enzymes were 
able to convert Dp 3-glc in Mv 3-glc with similar efficiency 
although in small proportion. 

Discussion 

Combining molecular physiology and genetics to dissect 
complex traits in plants 

Major steps forward in the study of plant anthocyanin 
synthesis were made through molecular physiology and 
forward genetics experiments [7]. However, QTL and 
association mapping techniques have also proven to be 
relevant tools in the dissection of complex physiological 
traits, in particular for anthocyanin metabolism [35]. Since 
their large scale application in maize, only a few attempts 
to clone and validate the polymorphisms of the genes 
underlying variation have been successful [36]. In this 
work, we propose different approaches to draw a global 
picture of the regulation of anthocyanin methylation in 
grape. As a starting point, we performed a rational sam- 
pling of highly diversified phenotypes present in the Vassal 
germplasm repository (INRA, Domaine de Vassal, France) 
by using the M-method [37,38] . This method ensures the 
optimal capture of phenotypic diversity suitable to detect 
contrasted genetic effects. The use of QTL mapping and 
the colocalization of a VvAOMr cluster on chromosome 
1, provided independent evidence supporting that the 
VvAOMT genes may play a key role in determining the 
level of methylated anthocyanins in grape. The availability 
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Figure 5 Difference in catalytic efficiency between the PB-A0MT2 and PK-A0MT2 allelic variants. Specific catalytic efficiency in 
performing 3'-and 3'5'-0-metliylation was measured on either cyanidin 3-0-glucoside (Cy 3-glc) and delpliinidin 3-0-glucoside (Dp 3-glc). Box 
plots indicate 5, 25, 50, 75 and 95 percentile, respectively. * indicates significant ANOVA test (P-value < 0.05). 



of broad genetic resources together with the use of grape 
genomic tools [31] helped refine our QTL detection 
through gene expression and association mapping studies, 
demonstrating the likely involvement of both c/5-regula- 
tion and newly characterized allelic variation in shaping a 
complex phenotype. 

The pleiotropic influence of VvMybA genes 

The QTL for MT level detected on LG2 (explaining 35% 
of the variance) co-localized with a VvMybA genes cluster, 
supporting that VvMybA genes have an important effect 
on anthocyanin methylation. The results found in previous 
work with the same experimental design [20] showed 
nonetheless that this region better explained the total 
anthocyanin content (64% of variance explained) than the 
level of methylated anthocyanins. This suggested that this 
QTL, cosegregating with VvMybA genes cluster, mostly 
had an effect on total anthocyanin accumulation and by 
consequence, an indirect effect on anthocyanin methyla- 
tion. The VvMybA genes control the transcription of key 
genes of the late steps of the anthocyanin biosynthesis: 
UFGT and AOMT [18]. Our study showed that VvMybAl 
has a significant effect on global VvAOMTs expression 



level. WZiFGr catalyzes the glycosyl transfer on C3 of the 
anthocyanidin aglycone [19,39]. WAOMr catalyzes the 
methylation of Dp 3-glc into both Pt 3-glc and Mv 3-glc, 
and that of Cy 3-glc into Pn 3-glc [11]. This pleiotropic 
effect of the VvMybA locus explains the relationship 
between the amount and the nature of the synthesized 
anthocyanins. This effect could also explain why black 
grape berries always contain large amounts of malvidin, 
the most methylated form of anthocyanidins [24-26]. 
UFGT and AOMT genes showed similar expression pro- 
files [18,40] and to date, polymorphisms in their sequence 
were never associated with variation in anthocyanin con- 
tent [9,11,13]. Thus, these genes were considered so far to 
be purely c/5-regulated and to have no regulatory control 
on the anthocyanin pathway. 

Determining the influence of VvAOMT genes 
independently of VvMybA genes 

We intended to identify regulation factors specific to the 
anthocyanin methylation level and independent of the 
total anthocyanin content. The MT variable and the 
total anthocyanin content were statistically independent, 
so the cosegregation of the QTLs for the two traits and 
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the correlation between the VvMybAl expression and 
MT were thus strictly due to biological effects. Nonethe- 
less, genes specifically inducing differences in anthocya- 
nin composition should either be independent of the 
expression of VvMybA or show trans -vegwldition or 
structural variation associated with anthocyanin methy- 
lation. In the search for such specific factors of the 
anthocyanin methylation process, the key was to find 
elements of variation that were significantly independent 
of the global anthocyanin VvMybA regulated chan- 
nel Recently, a similar example was shown in Ipomoea 
for a dihydroxyflavonol reductase (a gene also involved 
in the anthocyanin pathway) where duplicated genes 
escaped the control network for greater adaptation [41]. 
We might also suspect this to be the case of VvA0MT2 
and anthocyanin methylation. 

The VvAONlT genes cluster, the fate of a duplicated gene 
family 

Grape varieties accumulating only non-methylated 
anthocyanins have never been described [24-26] and 
were also not identified in our study on a very diverse 
sample, necessarily involving the persistency of functional 
AOMT in the grape genome. Two mechanisms can gen- 
erate such a pattern: a strong selective constraint on 
AOMT activity and a functional redundancy limiting the 
effect of deleterious mutations at WAOMriocus. We 
identified 3 candidate genes encoding OMTs located on 
chromosome 1 in the confidence interval of the QTL, 
potentially all coding for functional proteins, which sup- 
ports initial functional redundancy. In three genotypes 
with very contrasting MT phenotypes, only VvAOMTl 
and VvA0MT2 were expressed in mature grape berries, 
no expression of VvA0MT3 being detected. The total 
expression of VvAOMT {VvAOMTl and VvA0MT2 
together) showed a high correlation with the level of 
methylation, while the specific expression of VvA0MT2 
did not. We may thus hypothesize that VvAOMTl is cis- 
regulated by the VvMybA genes, as reported in Ageorges 
et al. [17] and Cutanda-Perez et al. [18], and broadly 
determines the level of methylated anthocyanins. Conver- 
sely, VvAOMTl expression was only marginally affected 
by VvMybA genes and showed structural variation, 
affecting anthocyanin methylation: the functional redun- 
dancy is decreasing as mutations accumulate. Although 
many other elements of variation may be involved in 
anthocyanin methylation, we may hypothesize a model 
as follows: VvAOMTl is the fundamental isogene for 
anthocyanin methylation under strict functional con- 
straint, while VvA0MT2 undergoing lesser functional 
constraint is responsible for fine and specific differences 
in the level of methylated anthocyanins (see proposed 
model Figure 6). 



The complex base of anthocyanin methylation in grape 

Previous works in petunia and grape suggested that dif- 
ferent OMTs were responsible for the 0-methylation of 
anthocyanins at either 3' or 5' positions [42,43]. In con- 
trast, Hugueney et al. [11] identified a bifunctional 
AOMT catalyzing 3' and 5' 0-methylation, suggesting 
that AOMT alone would be sufficient to perform both 
methylation steps in grape berries. In this study, we 
showed that the two proteins encoded by allelic variants 
of VvA0MT2 were also able to perform in vitro 3' or 
3',5' methylation on Cy 3-glc and Dp 3-glc but in addi- 
tion, we identified two SNPs associated with variation in 
the level of methylated anthocyanins. We showed that 
this polymorphism led to changes in substrate-specific 
catalytic efficiency, which only affects Dp 3-glc methyla- 
tion into Pt 3-glc. Differential 3' 0-methylation on trihy- 
droxylated anthocyanin substrate due to protein variants 
may explain the quantitative differences observed in 
planta, after removing the effect of VvAOMTl, Further- 
more, the comparison with the AOMTl and FAOMT 
characterized previously [11,12] showed that AOMT2 
variants characterized in this study tend to produce a 
lower proportion of Mv 3-glc, with a preferential cataly- 
sis of the 3' 0-methylation. Indeed, both AOMT2 var- 
iants synthesized the same low relative proportion of 
Mv 3-glc (8.1% and 5.3%, respectively) while in similar 
conditions, AOMTl was able to synthesize up to 34% of 
Mv 3-glc [11]. Finally, variation in catalytic properties 
together with potential differential expression is defi- 
nitely capable of creating the complex pattern of methy- 
lated anthocyanin variation we observe in different 
grapes cultivars. 

Conclusions 

The aim of this work was to identify and quantify the 
effect of molecular factors controlling the variation of 
anthocyanin methylation in grape by combining genetic 
mapping and functional study. Most previous studies 
regarding the accumulation of anthocyanin in grape have 
pointed to the main involvement of the VvMybA locus 
[16,18-21,44]. Our findings confirmed the effect of this 
locus, but we made an important step forward by reveal- 
ing the influence of novel factors other than Myb type 
genes involved in variation of anthocyanin composition. 
We quantitatively assessed the effect of the VvAOMT 
genes cluster on anthocyanin methylation. Furthermore, 
we identified new gene variants encoding proteins with 
slight differences in catalytic properties from other 
AOMTs described so far, showing allele-guided substrate 
specificity. This study improves our knowledge on how 
differences in methylated anthocyanin level may exist 
independently of total anthocyanin content. The relative 
independence of these two traits is possible because of 
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Figure 6 Synthetic scheme of the regulation of anthocyanin methylation in grape berry. The red arrows present the main influence of 
MybA protein c/'s-regulating AOMTl, the key enzyme involved in anthocyanin methylation. In addition, A0MT2, less influenced by MybA, trans- 
regulate anthocyanin methylation through allelic structural variation influencing enzymatic efficiency. Finally, A0MT3 expression was not detected 
in mature berries. 



the partial functional redundancy between the VvAOMT 
isogenes, among which VvA0MT2 "escaped" the control 
of VvMybA. Finally, with a relaxed functional constraint 
on VvA0MT2 while VvAOMTl is still ensuring the per- 
sistency of the activity, variation could arise, giving rise 
to novel phenotypes. 



Methods 

Plant materials 

The plant material consisted of three populations: one 
cross-derived mapping population for QTL mapping 
and two nested natural populations for the association 
study. The mapping population (SxG) was a Fl progeny 
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of 191 individuals from a reciprocal cross between clone 
73 of Syrah (S) and clone 516 of Grenache (G) repli- 
cated in two blocks (A and B) as previously described 
[20]. Fifty colorful genotypes from the agromorphologi- 
cal core-collection defined by Barnaud et al. [33] were 
selected to optimally represent the diversity in anthocya- 
nin methylation level. All genomic DNA analyses pre- 
sented in the study were realized on this collection of 
50 colorful cultivars. Within these 50 genotypes, we 
selected a subset of 32 individuals (listed in Additional 
file 2) designed through the M-method [37,38] to opti- 
mize the variability of the level of anthocyanins methyla- 
tion. All transcriptomic analyses and sequencing of 
cDNA performed in the study were realized on this 32 
subcore-coUection. Mature berries of each genotype 
were harvested 32 days after veraison. 

Traits analyzed 

Grape berries were harvested at the same level of matur- 
ity (20° Brix) in 2006 for the SxG progeny and in 2008 
for the colorful genotypes and the subcore-coUection. 
The anthocyanin composition of samples was analyzed as 
described in Fournand et al. [45]. The concentration of 
delphinidin (Dp), cyanidin (Cy), peonidin (Pn), petunidin 
(Pt) and malvidin (Mv) derivatives was expressed in mg 
of malvidin 3-0-glucoside equivalent per gram of fresh 
berry skin. Level of anthocyanin methylation (MT) was 
calculated as the ratio of methylated/non-methylated 
anthocyanins. The methylated/non-methylated ratio was 
log-transformed in order to unskew their distribution. In 
the SxG sample, the normality of the distribution was 
checked using the Shapiro- Wilks test for individuals con- 
taining anthocyanins (non-white cultivars). The total 
anthocyanin content was calculated as ln(l + x) where x 
is the sum of all anthocyanin compounds. The HPLC 
data and the MT variable for the subcore-collection are 
presented in Additional file 2. 

DNA extraction, genotyping and sequencing 

Two square centimeters (80-100 mg) discs of fresh young 
leaves were collected for each genotype. DNA was 
extracted using a Qiagen DNA Plant Mini Kit (QIAGEN 
S.A., Courtaboeuf, France) with slight modifications as 
described in Adam-Blondon et al. (2004). The 191 SxG 
offsprings have been previously genotyped for 97 SSR 
markers to cover the genome with minimal inter-marker 
space (maximum of 10 cM) with the most informative 
polymorphisms (priority given to 1:1:1:1 segregating mar- 
kers). An additional marker corresponding to a 200 bp 
insertion/deletion polymorphism (InDel) in the VvAOMTI 
locus promoter was genotyped. PGR conditions were 95°C 
for 7 min, 35 cycles of 95°C (30 s), 58°C (30 s), 72°C 
(1 min) and a final cycle of 72°C for 5 min. Amplified frag- 
ments were run on a 1% agarose gel, stained with 



ethidium-bromide and photographed under UV light. 
Amplification primers were designed using the Primer 3 
software and are listed in Additional file 10. PGR frag- 
ments were amplified, sequenced, and analyzed as 
described by Le Cunff et al. [46]. 

RNA extraction, sequencing and quantitative real-time 
PGR 

For transcriptomic assays, berries were sorted by density 
on salt gradient between 100 and 120 g/L to ensure homo- 
genous maturity stage within and across all samples and 
immediately frozen in liquid nitrogen. RNA was extracted 
from 200 mg of berry powder using RNeasy Plant Mini 
Kit (QIAGEN S.A., Courtaboeuf, France) following the 
manufacturer s protocol and then quantified with Ribo- 
green (Molecular probes, Leiden, The Netherlands). 
Reverse transcription was performed in triplicate from 
each sample from 500 ng of purified RNA as described in 
Fournier-Level et al. [20]. Amplification primers were 
designed using the Primer 3 software and are listed in 
Additional file 10. PGR fragments were amplified, 
sequenced and analyzed as described in Le Gunff et al. 
[46]. The expression of VvAOMT and VvA0MT2 was 
measured for the 32 cultivars showing contrasted methy- 
lated anthocyanin berry content. PGR was performed as 
described in Fernandez et al. [47] by comparing the cycle 
threshold (GT) of the target gene with that of VvEF-1 
alpha gene used as an internal standard. 

Framework genetic maps and QTL detection 

Framework genetic maps were constructed according to 
the pseudo-testcross strategy [48] using Garthagene 
0.999R software [49] as described in Fournier-Level et 
al. [20] with an additional marker added corresponding 
to a 200 bp InDel in VvA0MT2 of chromosome 1 seg- 
regating in the SxG cross. QTL detection was performed 
on both parents and consensus maps. Composite inter- 
val mapping (GIM) on MT was performed separately on 
the A and B blocks for the consensus map using 
MapQTL 4.0 [50], and both blocks together for parental 
maps using the multi-trait option QTL Cartographer 
[51]. For GIM on the parental maps, we used the For- 
ward & Backward regression method for cofactor selec- 
tion with 0.1 as the in-and-out threshold for the P- value 
of the partial F-test. Then, a genome scan was per- 
formed with a maximum of five cofactors within a 10 
cM window. LOD thresholds corresponding to an 
experimentwise error rate of 5% were then determined 
through 1000 permutations. For GIM on the consensus 
map, we determined the simple interval mapping LOD 
thresholds through 1000 permutations with a genome- 
wide error rate of 5%. 99% and 90% confidence intervals 
were determined through 1000 bootstrap resampling 
using the R/qtl package. 
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Statistical tests 

Statistical tests were performed using R software (The R 
project, USA) following the Shapiro- Wilks method for 
normality tests and Pearson's for correlation tests. Candi- 
date genes were selected on the grape reference genome 
sequence according to (1) their presence in the QTL inter- 
vals and (2) their homology with genes involved in antho- 
cyanin metabolism. Association tests were carried out 
using R following the Efficient Mixed Model Association 
method [34] with kinship matrix estimated using 20 SSR 
markers scattered along the 19 grape linkage groups [52]. 
Association tests were performed on each gene in order to 
determine the highly associated polymorphisms based on 
P-values and adjusted R ^ of a mixed linear model using 
the kinship matrix as the variance-covariance matrix to 
control for the genetic background "structure" effect. 

Cloning of VvA0MT2 and characterization of recombinant 
VvA0MT2 

The full-length cDNA of VvA0MT2 was amplified from 
cDNA of mature berries from two cultivars, Petit Bouschet 
X Aramon n°4 (PB) and Papadiko (PK), with high-fidelity 
Taq polymerase (Advantage-HF 2 PGR kit; Clonetech) 
using the upstream primer 5'-GTGGATCCTCTTCAAC- 
CATGTCCAGCTCAAGTCA-3' and the downstream pri- 
mer 5'-TTGCGGCCGCATAGATTTAGGCTAATAGA 
GGCGC-3'. The corresponding amplified cDNAs were 
named PK-AOMT2 (cv. Papadiko) and PB-OMT2 (cv. 
Petit Bouschet x Aramon n°4). Both amplified cDNAs for 
VvA0MT2 were cloned into the pGEM-T easy vector 
(Promega) and the resulting plasmids were sequenced to 
verify that no mutation had been introduced. The recom- 
binant plasmids were then digested by the BamHl/Notl 
enzymes, and the resulting restriction fragments were 
ligated into the pGEX-4 T-2 recombinant vector (GE 
Healthcare) digested by BamHl/Notl. The AOMT2 pro- 
teins were expressed as GST fusion proteins. Recombinant 
proteins were purified using glutathion Sepharose affinity 
resin following the pGEX instruction's manual. The char- 
acterization of the recombinant proteins was performed 
after the cleavage of the GST moiety using a site- specific 
protease whose recognition sequence was located immedi- 
ately upstream from the multiple cloning sites on the 
pGEX plasmid and then quantified by SDS-PAGE (Addi- 
tional file 11). The activity of both PB-AOMT2 and PK- 
AOMT2 enzymes was then tested in vitro using anthocya- 
nins as substrate in the presence of S-adenosyl-L-methio- 
nine (SAM) as described in Hugueney et al. [11]. The 
enzymatic reactions were incubated at 30°C for 1 h 30. For 
AOMT2 characterization, three independent recombinant 
protein preparations were tested, and four enzymatic 
assays were performed with each preparation. Reaction 
products were analyzed by HPLC-DAD as previously 
described [45]. Anthocyanins were identified according to 



their UV-visible absorption spectra and retention time; 
quantifications were based on peak areas at 520 nm using 
external calibration with malvidin 3-glucoside. 

Accession numbers 

Sequence data from this article have been deposited 
with the EMBL/GenBank data libraries under accession 
number # HQ702997 (AOMT2 from Petit Bouschet x 
Aramon n°4). The sequence Papadiko differs from Petit 
Bouschet X Aramon n°4 in two SNPs (T127 to C and 
A388 to T). 

Additional material 



Additional file 1: Summary statistics describing the anthocyanin 
profiles in the SxG progeny. The anthocyanin content is expressed in 
mg of aglycone anthocyanin per g of fresh berry skin. Key: Min, 
minimum; Qu, quantile; Max, maximum. 

Additional file 2: Anthocyanin content and expression profile of 
VvMybA, VvAOMT and VvA0IVIT2 in the Core-coll 32 cv. 

Anthocyanins content is expressed in mg of anthocyanin per g of fresh 
skin and expression profile for VvMybA, VvAOMT and VvA0MT2 
expressed as-2ACt. 

Additional file 3: Comparison of the MT variation between the two 
alleles of the SxG progeny measured in the block B. The Syrah 
parent is homozygote for the "strong" allele of VvOMT2 while the 
Grenache parent is heterozygous with the "weak" and"strong" 
combination of alleles. 

Additional file 4: List of the unigene from the NCBI Vitis EST 
database present in 1.05 Mbp surrounding the position of 
maximum LOD for the QTL on the LG1. 

Additional file 5: Comparison of the A0MT1, A0MT2 and A0MT3 
amino acid sequences. Amino acid sequences were aligned using 
CLUSTAL W. Residues identical to the VvAOMTl sequence were shaded. 

Additional file 6: Specific amplification of AOMT, A0IVIT2 and 
A0I\/IT3 on both (A) genomic DNA and (B) coding DNA. (A) 

Amplification on genomic DNA (cv. Syrah). (B) Amplification on cDNA of 
cv. Joubertin (1; 2; 3), cv. Molinara (4; 5; 6) and cv. Petit Bouschet (7; 8; 9) 
mature berries of AOMT (1; 4; 7), A0MT2 (2; 5; 8) and A0MT3 (3; 6; 9). 

Additional file 7: Linkage disequilibrium between polymorphisms in 
the VvAOMTl and VvA0IVIT2 genes for the collection of 50 cv. R2 

are reported on the lower diagonal and P-val from chi-square test for 
independence are reported on the upper diagonal. 

Additional file 8: Quantitative real-time PCR expression profiling of 
AOMT (black) and A0MT2 (grey) in hairy roots (control and 
VlmybA1-2 transformed). All data are the mean (± SD) of three 
replicates on two (control) and four (VlmybAl-2 transformed) 
independent lines. Expression values have been normalized with 
VvEFl alpha and expressed as relative abundance. 

Additional file 9: Polymorphisms in the VvA0MT2 cDNA. The name 
of the SNPs corresponds to the nature of the mutation in lUB code 
followed by the position on the coding sequence. Non-synonymous 
mutations are presented in bold with a red arrow, mutations in the 
3'UTR non-coding sequence are presented with dashed arrows and not 
reported at the right scale. 

Additional file 10: PCR primers used for real-time PCR and expected 
size for amplified fragments. Key; gDNA (genomic DNA); cDNA (coding 
DNA). 

Additional file 11: SDS-PAGE of the recombinant AOMTs obtained 
from cv. Papadiko (PK-AOMT) and cv. Petit Bouschet x Aramon n°4 (PB- 
AOMT). Lanes: (1) protein size markers; (2) purified PK-AOMT after GST 
cleavage. (3) purified PB-AOMT after GST cleavage. The molecular 
weights of the markers are indicated in kDa. 
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